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Abstract 


Using the efficient discrete- ordinate method, we present an analytical solution for 
radiative transfer in the coupled atmosphere-ocean system with rough air-water interface. 
The theoretical formulations of the radiative transfer equation and solution are described. 
The effects of surface roughness on radiation field in the atmosphere and ocean are 
studied and compared with measurements. The results show that ocean surface roughness 
has significant effects on the upwelling radiation in the atmosphere and the downwelling 
radiation in the ocean. As wind speed increases, the angular domain of sunglint broadens, 
the surface albedo decreases, and the transmission to ocean increases. The downward 
radiance field in the upper ocean is highly anisotropic, but this anisotropy decreases 
rapidly as surface wind increases and as depth in ocean increases. The effects of surface 
roughness on radiation also depend greatly on both wavelength and angle of incidence 
(i.e., solar elevation); these effects are significantly smaller throughout the spectrum at 
high sun. The model-observation discrepancies may indicate that the Cox-Munk surface 
roughness model is not sufficient for high wind conditions. 


OCIS codes: 010.1290, 010,4450, 030.5620 
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1. Introduction 


Cox and Munk 1 described the statistical charactersitics of reflection by wind-blown ocean 
waves by modeling the sea surface as a collection of individual mirror facets. They 
presented the probability distribution for the slopes of surface facets as a wind-speed 
dependent Guassian function. Based on this Cox and Munk formulation, several 
researchers incorporated the ocean surface roughness in their radiative transfer models ' . 
Most of these models used the ray tracing method or the Monte Carlo technique to treat 
the surface roughness. Monte Carlo approach consists of using probabilistic concepts and 
has the advantage for geometries other than the plan-parallel. Implementation of the 
statistical surface roughness by Monte Carlo method is relatively straightforward. The 
discrete-ordinate technique, on the other hand, can be more computationally efficient and 
accurate, because it solves the radiative transfer equation analytically without the 
enormous statistical sample required to close a Monte Carlo solution and without the 
statistical fluctuation error. However, due to its analytical nature, implementation of the 
surface roughness in a discrete-ordinate radiative transfer model is more complicated; a 
rigorous solution involves an additional parameter that results in a different analytical 
solution from the flat surface case. 

To extend applications of the DIScrete-Ordinate Radiative Transfer (DISORT) code for 
systems including two media (atmosphere and ocean, atmospthere and ice, etc.), Jin and 
Stamnes 8 developed a Coupled DISORT (CDISORT). The CDISORT code accounts for 
change in the refractive index change at the boundary of the two media. For radiative 
transfer in such a coupled system, CDISORT treats the ocean or ice same as atmospheric 
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layers but with different optical properties, particularly, different refractive indices. 
However, the interface between two strata with different refractive indices was 
considered as flat. This flat surface assumption limits the applications of CDISORT; the 
afore-mentioned wind-blown ocean surface is hardly flat. In addition to affecting 
reflection, the surface roughness itself significantly affects the directional character of the 
beam transmitted beneath the air- water interface. Gjerstal and co-workers 9 proposed an 
ad hoc method to consider the surface roughness in the discrete-ordinate method. They 
mimic the irradiances from a Monte Carlo model by adjusting the refractive index in 
CDISORT . This method has a number of limitations, for example, it calculates 
irradiances in the ocean only. In this paper, we present a more consistent and widely 
applicable solution of the discrete ordinate radiative transfer problem in the coupled 
atmosphere-ocean system with rough surface. 

2. Equation and Solution of Radiative Transfer 

In order to incorporate the ocean surface roughness into the radiative transfer equation 
and obtain an analytical solution through the discrete- ordinate method, we need to make 
the following assumptions: 

> The rough surface can be resolved as a series of small planar facets and, the 
orientations (slopes) of these facets follow a certain statistical distribution, for 
example, the Gaussian distribution described by Cox and Munk 1 . 

> The dimensions of the elemental facets and surface undulations are large 
compared with the wavelength of light, so geometric optics can be applied to 
calculate the reflection and refraction at the surface. 
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> The optical depth of either ocean or atmosphere is independent of the surface 
roughness or horizontal position because statistically, there is no difference 
between any two points on surface. 

Under these assumptions, the time-averaged radiative effects at any two points on surface 
are the same, and a patch of surface area at an instant in time in which every possible 
slope occurs can represent the surface as a whole. Therefore, radiative transfer in a 
coupled system with horizontally homogeneous atmosphere and ocean and with rough 
ocean surface is still in the one-dimensional category, as long as the calculated radiation 
is considered to be time-averaged (statistical averaged) for a point or to be relevant to 
surface area larger than the patch afore-mentioned for a instant. 

o 

Jin and Stamnes (hereafter refers as JS94) presented in detail the solution for coupled 
(i.e., air-sea) radiative transfer by the discrete-ordinate method for the flat ocean case. 
Here we will follow the same conventions defined in JS94 to describe the discrete- 
ordinate radiative transfer equation and solution for the rough surface case. Because the 
formulations have a lot in common between these two cases, we will omit most common 
derivations and emphasize the differences here. 

For the flat surface case, the ocean was divided into a totally reflecting angular domain 
where upwelling photons cannot return directly to the atmosphere, and a refracting 
domain where upwelling photons can pass the interface directly to the atmosphere. 
However, once the surface roughness is introduced, there are no such distinct angular 
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domains. Because of the possible ranges of angles for both the incident photons and the 
surface interfaces, photons in the atmosphere may pass an interface directed to any angle 
downward; and vice versa for photons from ocean to atmosphere. This difference results 
in different radiative transfer solutions. 

The radiative transfer equation to be solved for a plane-parallel medium with one 
dimension can be written as 

li ^ = J(7’ </>) ~ C ^~ Pi*, f ) I( f, M' , </>' W + Q(t, </>) (1) 

where I(t,h,</> ) is the radiance at vertical optical depth x (measured downward from the 
upper boundary) and in direction (p, (f>) ; pt is the cosine of the zenith angle (positive for 
upward directions); 0 is the azimuth angle; co is the single scattering 
albedo; p(T,jU,(/),jU' ,(/)') and Q(t,jU,0) are the phase function and source function 
respectively. We only consider the solar radiation (i.e., not terrestrially emitted thermal 
infrared or microwave). The solar beam source, , is different from the case of a 

flat ocean. In the flat surface case, part of the down welling solar beam is reflected 
specularly back to the atmosphere, and the rest is refracted into the ocean at an angle 
which depends on the refractive index; this results in two terms (downwelling and 
reflected) in the solar source function for the atmosphere (Equation 3 in JS94) and one 
refractive index dependent term for the ocean (Equation 4 in JS94). However, for the case 
of rough ocean, the solar beam is diffused to various directions when it hits the surface. 
Therefore, there is no beam source term in the ocean and only one expression in the 
atmosphere for the rough ocean case, which is 
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,,, , om i)/ ^ p( r u 4 ) _y^ ^ , cxp( _ T f fl ^ } r < r, 

S(r, //.« = { Q , /4 * r>r _ (2) 

where x a is the total optical depth of the atmosphere, p () is the cosine of the solar zenith 
angle, (f>o is the solar azimuth angle, and Fo is the solar-beam intensity at the top of the 
atmosphere. 


Expanding the radiance I(t,jU,(/>) into a Fourier cosine series of 2N and the phase 


function p(T,jU,(/),jU' ,</>') into a series of 2N Legendre polynomials, the discrete-ordinate 
method converts Equation (1) into a system of azimuthally independent, coupled 
differential equations for each of the Fourier components. Detailed derivations of these 
equations were given in JS94 and will not be repeated here. Following the same 
procedure, the equations for each azimuth radiance component (here we omit the index 
denoting the order of Fourier series) can be derived, which are in the atmosphere 


Pi 


dr 


= I(t, jU“)~ 2 co) D(t, h* , jU) )/ ( r, ju)) + X 0 (T , jU) ) exp(- r/// 0 ) , 


j=-Ni 
i* o 


(3a) 


/ = ±l,±2,---,±2V i; r < r a 


and in the ocean 




dl{r,p°) 

dr 


= ^CQ°D(r, /i° ,fi°) I 


j=-N 2 

j* 0 


(3b) 


i = ±l,±2,---,±N 2 ; t > T a . 

Equations (3a) and (3b) are analogous to the Equations (7) and (8) in JS94, but with 
different source terms. Here 2Ni and 2 Nt are the numbers of quadrature points (i.e., 
stream numbers) applied in the atmosphere and ocean, respectively. D{r,jU j . // , ) and 


X 0 (t, fli) were also defined in JS94. The ( //," , of) ) and ( //," . a ) ) are quadrature points and 
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weights for the atmosphere and ocean, respectively, with fi_. = — // ■ and = co i . They 


have the following relationships: 


Mi 


^-[l-l//;) 2 ]/^) 2 , 1 = 1,2,-,^, 


(4) 


co 




ty. 


i = 1,2, A/',. 


(5) 


Here n a and n 0 represent the refractive indices of air and water, respectively. Following 
the same procedure as in JS94, the solutions for Equations (3a) and (3b) can be obtained 


as 


/( r, //, fl ) = X [c_, G_ . (ji“ ) exp(^; r) + C, G ; (//,« ) exp(-k; t)\ +Z 0 (<u“ ) exp(- t / h 0 ) , 


7=1 


(6a) 


/ = ±1,±2,---,±A^,; r < r fl 


= [C ,G_, (//" ) exp(A:;T) + C y G (//,° ) exp(-£; r)l 

i = ±1,±2, ■■•,±iV 2 ; r>r fl . 


7=1 


(6b) 


with constants C +/ . k } and G f are eigenvalues and eigenvectors, respectively, 
determined by solving an algebraic eigenvalue problem as described in Stamnes et al. 10 . 
Z 0 ( // " ) is defined and obtained by solving the Equation (12b) in JS94. The solutions 
represented by (6a) and (6b) are seemingly simpler than those in JS94 for the flat ocean 
case. However, the solution is not complete yet, because constants C ±J in (6a) and (6b) 

are still unknowns which differ from layer to layer in the atmosphere and the ocean (for 
simplicity, we omitted the index denoting layers here). These constants will be 
determined by boundary and interface conditions for radiances (intensities). 
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The conditions for the top and bottom boundaries, for the interfaces among atmospheric 
layers, and for the interfaces among oceanic layers are same as those for the flat ocean 
case, which were given by Equations (16a), (16b), (16f) and (16g) in JS94. However, the 
continuity conditions for radiances at the interface between the atmosphere and ocean are 
very different from those for the flat ocean case. If we denote t~ as the optical depth just 


above the ocean surface and T + a as that just below the surface, these conditions for the 
rough surface case can be expressed as 


rvi 2 

) = X f(T ‘‘ -Vj ’ n ° ) + X I(T+a ' M ° j )r( ^ ’ V°j ’ n a K ) 

;■= 1 7=1 


+ ~Mo F o exp(-T fl /jU Q ) R(fi? -JU 0 ,n a /n a ), (7) 

n 

i = 1,2 iVj 


V, n 2 


7=1 7=1 

+ -/Jo F o exp (-T a /jU 0 ) T (-A,° -Mo . n o / n a X (8) 

71 

i = l,2,-~,N 2 

in which reflection R and transmission T matrices appear without indices denoting the 
Fourier order. Equations (7) and (8) show that the emerging radiance at any direction at 
the air- water interface depends on incidences from all directions from both the 
atmosphere and ocean for the rough ocean case. This contrasts with the simple one to one 
correspondence (pairing of each //, a and //,°) as presented by equations (16c- 16e) in JS94 
for the radiances across the air-water interface of a flat ocean. One asset of the rough 
ocean case is the term accounting for the diffusion of the solar beam (the last term in (7) 
and (8)); it makes a simpler formulation of the particular solution (Equations (6a) and 
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(6b)) than that for the flat ocean case in JS94. The reflection and transmission matrices 


are calculated as 


r\ ^ 

R(jU,ju',n) = </>-</>', n,ju n )cosm(<p-</>')d(</)-</>') 


T(ju,ju',n)= 2 n , jU n ) cos m(<p - </>')d - </)') 

JO 


In Jo 


8 = ( *’ 
U m 0 l Q 


if m - 0 
otherwise 


(9) 


( 10 ) 


Here R(ju,ju',<p-<p',n,ju n ) and T (ju,ju\ <f>- </>',n,jU n ) are the reflection and transmission 
functions 7 ’ 11 , respectively. (p\ ()’) and (p, (|)) represent the incident and exit light 
directions, respectively. Note that in these functions, n is the relative refractive index, 
which equals n 0 /n a if the incident light is from the air (n a /n 0 if the incidence is from 

the ocean). p n is the cosine of the normal to the surface facet, and it follows a probability 
distribution which defines the surface roughness. p n is related to p and p' as 

// + //' 


Mn = 




( 11 ) 


a = )l — // 2 )l- ju' 2 cos(tp-tp') 


( 12 ) 


Substituting (6a) and (6b) into the boundary and interface conditions [i.e., the Equations 
(16a), (16b), (16f) and (16g) in JS94 and the Equations (7) and (8) here], we obtain a 
system of linear algebraic equations for the unknown coefficients Cj. The method to solve 
the equations and obtain the unknown coefficients was described in Stamnes et al. 10 and 
is not repeated here. The implementation of these solutions into the CDISORT code is 
not trivial, however. 
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3. Examples of Model Simulations 

3.1. Brief Description of Model 

The CDISORT just described has been used as the radiative transfer solver by our 
Coupled Ocean-Atmosphere Radiative Transfer (COART) model 1 ' 14 ( http://www- 
cave.larc.nasa.gov/cave/) . CDISORT accounts for the change in refractive index at the 

o 

air-sea interface and now includes the interface roughness into the analytic solution of 
the radiative transfer equation. Hence COART considers the atmosphere and ocean as 
one system and treats the ocean strata just as additional “atmospheric” layers with 
different optical properties. CO ART models the absorption and scattering processes in 
atmosphere and ocean explicitly. These include the scattering and absorption by 
molecules, aerosols and clouds in the atmosphere, and by liquid water molecules, 
dissolved and particulate matter in the ocean. 

COART calculates radiances and irradiances at any level of the atmosphere and ocean in 
both narrowband (spectral) and broadband. For the narrowband scheme, users can specify 
both the band (wavelength) limits and computational resolution arbitrarily. In this 
scheme, COART employs the LOWTRAN 7 band model (spectral resolution of 20 cm 1 ) 
and molecular absorption database for the atmosphere. This corresponds to a wavelength 
resolution of about 0.5 nm at 500 nm and 8 nm at 2000 nm. For efficient broadband 
calculations of radiance and irradiance, COART divides the solar spectrum (0.20-4. 0pm) 
into 26 fixed wavelength intervals; in each spectral interval, the k-distribution technique 
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parameterizes molecular absorption in the atmosphere using the HITRAN 2000 
database 21 . 

The most prominent effects of ocean surface roughness on solar radiation are on 
upwelling fields in the atmosphere and downwelling fields in the ocean. Its effects on the 
downwelling radiation in the atmosphere and the upwelling radiation in the ocean are 
significantly smaller. While COART can simulate a variety of quantities, including the 
water-leaving radiance, we show here mainly the types of calculations that pertain to 
surface roughness. 

3.2. Effects of Surface Roughness on Radiance 
We use the Cox-Munk 1 probability distribution of surface slopes for a given wind speed 
to model the roughness of the ocean surface. The slope distribution is 

P{ tan0„) = exp( ) (13) 

where 9 n is the tilt angle of the surface facet normal. The mean slope distribution width 
G is related to the wind speed U (m/s) as 

< 7 2 =0.003 + 0.0051217 (14) 

The CO ART calculations in Figure 1, which use a McClatchey Midlatitude Summer 
atmosphere 15 with marine aerosol optical depth of 0.1 (at 500 nm) and Case 1 water for 

1 Z O 

ocean with chlorophyll concentration of O.lmg/m , span the upwelling radiance 
distribution at the top of atmosphere (TO A) and the downwelling radiance distribution at 
10m and 50m in the ocean for three different wind speeds (3, 6 and 10 m/s). Figure 1 uses 
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polar coordinates, with view zenith angle (0) on the radial axis and relative azimuth angle 
(()>) as the azimuthal coordinate. To facilitate comparisons of different wind speeds, 
wavelengths and levels, the radiance in Figure 1 is normalized by the upwelling or 
downwelling irradiance (E) at the same level to the Anisotropic Radiance Function ( ARF ) 
as 

ARF{0,(j)) = n1 ^ 6 ^ (15) 

E 

where E is the upwelling irradiance if 7(0,( |>) is the upwelling radiance (0 is positive in 
this case in Figure 1) and otherwise, E is the downwelling irradiance (0 is negative in this 
case). ARF here is in fact the ratio of the actual radiance, 7(0, (f>), and the imagined 
isotropic radiance, El n, with the same irradiance. Therefore, the gradient in ARF 
represents the departure of radiance field from the isotropic case (ARF=1.0). 

The solar zenith angle in Figure 1 is 40 degrees. Because the slope distribution in 
Equation (13) is independent of the wind direction, the ARF (and the radiance field itself) 
of Figure 1 is symmetric with the principal plane (the vertical plane containing the sun, 
the surface target and the nadir). So only ARF for azimuth from 0° to 180° is presented. 
The entire prinicipal plane is covered by the horizontal axis of each panel in Figure 1 and 
the sun (observer) is on the left (right). The top 3 rows of Figure 1 show upwelling ARF 
for the broadband shortwave (0.20-4. 0pm), 531 nm (the central wavelength of MODIS 
channel 11), and 865 nm (for MODIS channel 16), respectively. The bottom 4 rows 
show downwelling ARF, for broadband and 531 nm, at 10 m and 50 m below the ocean 
surface (865 nm is not shown because of strong absorption by liquid water). The 3 
columns of Figure 1 cover wind speeds of 3 m/s, 6 m/s, and 10 m/s. The hot spot in each 
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panel represents the specular reflection (the sunglint in atmosphere) or transmission (in 
ocean) of the solar beam at the rough surface. The sunglint is conspicuous at the right of 
each TOA panel (top 3 rows), but the sunglint region widens as wind speed increases and 
is much more prominent for the 865 nm because of less atmospheric scattering. And at 
ocean depths of 10 m and 50 m, the sharp peaks in downwelling ARF show that the 
anisotropy for downwelling is much more prominent than for upwelling at the TOA; this 
occurs mainly because of the lager refractive index of water than air (i.e., the refraction 
effect). Glancing at the middle rows 4-5 of Figure 1 for downwelling ARF at 10 m, the 
increase in wind speed (surface roughness) appears to have only minor impacts on the 
beam refracted from the interface; when we concentrate on just the forward scattering 
component of the principle plane in the next paragraph (and Figure 2), the impact of 
surface roughness on downwelling will be more apparent. As depth is increased to 50 m 
(bottom rows 6-7 in Figure 1), the diffusion by ocean optics becomes more important. 

For the same atmospheric and ocean inputs, Figure 2 further shows the radiance 
distribution at 531 nm only; and just around the hot spots in Figure 1 in the principle 
plane, where radiance varies most sharply and wind effect is most glaring. Figure 2 
highlights the different impacts of wind on the radiance fields at various levels in the 
atmosphere and ocean. While the color scale of row 5 in Figure 1 revealed virtually no 
effect of surface roughness on downwelling radiance broadly over the hemisphere, rows 
2-3 of Figure 2 show that wind indeed has an impact on the radiance distribution around 
the forward scattering direction in ocean. Figure 2 further delineates how the wind effect 
diminishes, and the radiance anisotropy rapidly decreases, as depth in ocean increases. 
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Figure 3 shows a model-observation comparison of the shortwave radiances at the TO A. 
The measurement data were from NASA’s Cloud and Earth Radiation Energy System 
(CERES) instrument 17 during a special field program at the CERES Ocean Validation 
Experiment (COVE) site . CERES was programmed to a special mode for intense 
observation at COVE and only the measurements in those clear days during the field 
experiment are presented here. In this experiment, comprehensive measurements on a 
variety of physical and optical properties of the atmosphere, surface and ocean were also 
available for the model input here. The horizontal coordinate in Figure 3 is the sunglint 
angle, defined as the angle between the view direction and the specular solar reflecting 
direction for an imagined flat surface. The nine numbers in the lower portion of Figure 3 
are the mean model-observation biases for the nine glint-angle intervals (10 degrees each) 
from 0 to 90, respectively. Though the aerosol loadings and the surface and ocean 
properties were different for different days, the model and observation agrees fairly well 
away from the sunglint region. The difference is somewhat larger near the sunglint center 
(smaller glint angles), probably due to the error in the surface roughness treatment in the 
calculations, for example, the uncertainties in the Cox-Munk model. The SZA is around 
20 degrees when CERES made the measurements, and so a large glint-angle (larger than 
75 in Figure 3) also represents a large view zenith angle, where the view path is longer 
and surface footprint is larger, and therefore the possible horizontal variations of aerosol 
and surface have larger effects than at a small view angle. This might be responsible for 
the increased biases in the large angle regime. 


13 



3.3. Effects of Surface Roughness on Irradiance and Albedo 


The effects of ocean surface roughness on irradiances are shown in Figure 4, which has 
upwelling irradiances in the atmosphere (linear scale) and downwelling irradiances (log 
scale) in the ocean for 531 nm, 865 nm and the broadband shortwave (the three columns) 
at four different levels (the four rows). The model inputs for the atmosphere and ocean 
are identical in Figures 1 and 4. In each panel of Figure 4, the irradiances for different 
wind speeds are presented as a function of the cosine of SZA. Results for a flat ocean 
case (wind = 0 m/s) are plotted as the solid lines in each panel, and thus the difference of 
irradiances between a rough ocean case (represented by an non-zero wind) and the flat 
ocean case represents the surface roughness effect quantitatively. Figure 4 shows that the 
effect of surface roughness is smaller for high sun than for low sun; and the upwelling 
irradiance just above the surface (row 2) is the field with the most action. Note that for 
upwelling irradiance at low sun, roughness has a larger effect at 865 nm (where much of 
the downwelling is direct, striking the surface at glancing angle, thereby obtaining much 
Fresnel reflection) than at 531 nm (where the downwelling is more diffuse and has a 
component that is closer to normal). However, we find the opposite for high sun: 
roughness has a larger effect on above-surface upwelling at 531 nm than at 865 nm. The 
turning point is at approximately SZA of 60 (cosSZA of 0.5), where the direct and diffuse 
reflectances are similar and this angle can be considered as an “effective” angle for 
diffuse radiation. For each rough surface case, the variation with respect to the flat 
surface case (the solid line) in the downwelling irradiance just below the surface (row 3) 
is equivalent to the variation in the upwelling irradiance just below the surface (row 2) 
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but with opposite sign. However, because the downwelling irradiance just below the 
surface is much larger than the up welling irradiance just above the surface, the relative 
variation in the downwelling irradiance in ocean from the flat ocean case to a rough 
surface case is smaller and is less obvious than in the upwelling in atmosphere in Figure 
4, especially for high sun conditions. At depth of 10 m (row 4), the effects of internal 
ocean optics on irradiance outweigh the effects of surface roughness and the irradiance at 
865 nm is none due to strong water absorption. 

Because the downwelling irradiance in the atmosphere has little dependence on the 
surface condition of an ice-free ocean, the large effect of surface roughness on upwelling 
energy (top half of Figure 4) will have a signal in the surface albedo. The left panel of 
Figure 5 shows the MFRSR measured (670 nm) (the dots) and modeled (the solid lines) 
surface albedo for three clear afternoons with quite different wind regimes (right panel) at 
COVE; the aerosol loadings were low. Aerosol optical properties used in model were 
measured from the same platform by NASA’s Aeronet Cimel instrument 19 . The Cimel 
sun photometer made periodic scans in the almacantur and in the solar principal plane; 
inversions of these data yielded aerosol phase functions and particle size distributions” . 
The wind data were from the NOAA meteorology station also at COVE. The ocean 
optical properties and chlorophyll concentration were also from in situ measurements 14 , 
but ocean optics has little effect on the total surface albedo at 670 nm. To remove the 
relative difference between the two surface-based MFRSR instruments and obtain 
accurate ocean albedos, the instruments subsequently used for the downwelling and 
upwelling spectral irradiance measurements were calibrated relative to each other in 
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advance, by observing the same target at the same time 14 ' 18 . Results in Figure 5 show the 
significant effects of wind on ocean surface albedo, especially for large SZA. The 
dependences of albedos on SZA and wind are consistent between model and 
measurements. 

When light is incident on the rough surface at a grazing direction, the photons are more 
likely to undergo multiple reflection among the surface wave facets. There is also a 
shadowing effect of one wave facet blocking rays from getting to another facet 
(occultation) . Figure 6 shows the effects of multiple reflection and shadowing among 
the surface wave facets on albedo simulation. In each panel, the solid line is the modeled 
albedo same as shown in Figure 5, with both shadowing and multi-reflection considered. 
The short dashed line is the calculation without shadowing but with multi-reflection, 
while the long dashed line represents the results without any multi-reflection but with 
shadowing considered. The error bar shows the range of measured albedo within 3 
degrees of SZA centered by the error bar. As expected, the calculated albedo is reduced 
after including the shadowing effect and is increased after including the multi-reflections. 
The effect of multi-reflection is larger on day 1 when the wind was high and the model- 
observation agreement is improved after the effects of shadowing and multi-reflection are 
taken into account. However, the effects of shadowing and multi-reflection are small for 
high sun (small SZA). 

Results above show that the surface roughness has largest effect on ocean albedo at low 
sun (large SZA). Both Figure 5 and Figure 6 show a larger model-observation 
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discrepancy for large SZA (higher than about 80) for the day with strongest wind (day 1). 
This may indicate a larger error tendency in the Cox-Munk surface roughness model in 
high wind conditions. An alternate distribution was recently produced by Ebuchi and 
Kizu~ , based on approximately 30 million satellite observations over five years. The 
Ebuchi-Kizu function has a narrower slope distribution, and less sensitivity to wind, than 
the Cox-Munk function. The calculated albedo (not shown) based on the Ebuchi-Kizu 
model is close to that based on the Cox-Munk model for low winds, but higher for high 
wind and large SZA. However, we can not conclude the superiority of either the surface 
roughness parameterization, based on limited comparison. Validation of this 
parameterization is not a subject here. We simply indicate that the surface slope 
distribution function has significant impact on the albedo calculation too. 

4. Conclusion 

We present an analytical approach for radiative transfer in a coupled atmosphere-ocean 
system having a rough surface between two media with differing indices of refraction. 
The discrete-ordinate technique is used in the formulation and solution. The solution is 
implemented in the radiative transfer code - CDISORT. Using CDISORT as the radiative 
transfer solver, a coupled ocean-atmosphere radiative transfer (COART) model is now 
available to calculate various radiances and irradiances at any altitude in the atmosphere 
and depth in the ocean. This model is demonstrated online at http://www- 
cave.larc.nasa.gov/cave/ . 
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Model simulations show that the ocean surface roughness has significant effects on the 
upwelling radiation in the atmosphere and the downwelling radiation in the ocean. As 
wind speed increases, the angular domain of sunglint broadens, the surface albedo 
decreases, and transmission through the air- water interface to the interior of the ocean 
increases. The transmitted radiance just below the ocean surface is highly anisotropic, 
but this anisotropy decreases rapidly as surface wind increases. Deeper below the 
surface, anisotropy decreases as the optical properties of the ocean interior eventually 
overcome the impact of surface roughness. The effects of surface roughness on radiation 
depend greatly on both wavelength and angle of incidence (i.e., solar elevation); these 
effects are significantly smaller throughout the spectrum at high sun. 

The models and observations agree fairly well on the effects of surface roughness. Some 
discrepancies may indicate that the original Cox-Munk surface roughness model is not 
sufficient for high wind conditions. 
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Figure caption: 

Figure 1. Model simulated upwelling radiance field at the top of atmosphere (TO A) and 
the downwelling radiance field at depths of 10 m and 50 m in the ocean for wind 
speeds of 3, 6 and 10 m/s and for three wavelength sets (broadband, narrowband at 
531 nm, and narrowband at 865 nm). The SZA is 40 degrees. 

Figure 2. Effects of surface roughness on radiance distributions at 531 nm in the 

components of the principal plane containing most of the reflected solar beam in the 
atmosphere (top row), and most of the refracted solar beam in the ocean (middle and 
bottom rows). 
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Figure 3. Comparison of modeled and measured broadband radiances as a function of 
sun-glint angle. The nine numbers are the mean model- observation biases for the nine 
glint-angle intervals (10 degrees each) from 0 to 90, respectively. 

Figure 4. Modeled irradiances versus cos(SZA) with upwelling irradiance in the 

atmosphere and downwelling irradiance in the ocean, for different wind speeds and 
different wavelengths. 

Figure 5. Effects of wind speed on ocean surface albedo at 670 nm. The left panel shows 
the modeled and measured surface albedo during three afternoons. The right panel 
shows the observed wind speed for each afterenoon. Different colors are for different 
days. 

Figure 6. Effects of multiple reflection and shadowing among surface wave facets on 
ocean albedo simulation. The three panels are for the three selected days as in Figure 
5. The long dashed line is the albedo computed without multi-reflection but with 
shadowing; the short dashed line is the albedo computed without shadowing but with 
multi-reflection. The error bar represents the measured albedo range within the 3 
degrees of SZA. 
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Figure 1 . Model simulated upwelling radiance field at the top of atmosphere (TOA) and 
the downwelling radiance field at depths of 10 m and 50 m in the ocean for wind 
speeds of 3, 6 and 10 m/s and for three wavelength sets (broadband, narrowband at 
531 nm, and narrowband at 865 nm). The SZA is 40 degrees. 
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Figure 2. Effects of surface roughness on radiance distributions at 531 nm in the 

components of the principal plane containing most of the reflected solar beam in the 
atmosphere (top row), and most of the refracted solar beam in the ocean (middle and 
bottom rows). 
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Figure 3. Comparison of modeled and measured broadband radiances as a function of 
sun-glint angle. The nine numbers are the mean model-observation biases for the nine 
glint-angle intervals (10 degrees each) from 0 to 90, respectively. 
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Figure 4. Modeled irradiances versus cos(SZA) with upwelling irradiance in the 

atmosphere and downwelling irradiance in the ocean, for different wind speeds and 
different wavelengths. 
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Figure 5. Effects of wind speed on ocean surface albedo at 670 nm. The left panel shows 
the modeled and measured surface albedo during three afternoons. The right panel 
shows the observed wind speed for each afterenoon. Different colors are for different 
days. 
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Figure 6. Effects of multiple reflection and shadowing among surface wave facets on 
ocean albedo simulation. The three panels are for the three selected days as in Figure 5. 
The long dashed line is the albedo computed without multi-reflection but with 
shadowing; the short dashed line is the albedo computed without shadowing but with 
multi-reflection. The error bar represents the measured albedo range within the 3 degrees 

ofSZA. 
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